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ABSTRACT 

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic 
Legacy Survey is to track galaxy evolution back to z « 8. Its two-tiered "wide and deep" strategy 
bridges significant gaps in existing near-infrared surveys. Here we report on z « 8 galaxy candidates 
selected as F105W-band dropouts in one of its deep fields, which covers 50.1 arcmin^ to 4 ks depth 
in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. 
Two of our candidates have J < 26.2 mag, and are > 1 mag brighter than any previously known 
F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity 
function of galaxies at z w 8, and show that the number density of such very bright objects is higher 
than expected from the previous Schechter luminosity function estimates at this redshift. Another 
two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first 
such individual detections at z w 8. Their derived stellar masses are on the order of a few xlO^ Mq, 
from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function 
at z « 8. The high number density of very luminous and very massive galaxies at z w 8, if real, could 
imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early 
time. 

Subject headings: cosmology: observations — galaxies: luminosity function, mass function 



> 
o 



^ Department of Physics & Astronomy, University of Missouri, 
Columbia, MO 65211, USA 

^ Department of Astronomy, The University of Texas at Austin, 
Austin, TX 78712, USA 

^ Department of Physics & Astronomy, Johns Hopkins Univer- 
sity, 3400 N. Charles Street, Baltimore, MD, USA, 21218 

Space Telescope Science Institute, 3700 San Martin Drive, Bal- 
timore, MD 21218, USA 

^ National Optical Astronomy Observatory, 950 N. Cherry Ave., 
Tucson, AZ 85719, USA 

^ Department of Physics and Astronomy, University of Pitts- 
burgh, 3941 O'Hara Street, Pittsburgh, PA 15260, USA 

Physics & Astronomy Department, Rutgers University, Piscat- 
away, NJ 08854, USA 

Astronomy Department, The University of Arizona, Tucson, 
AZ 85721, USA 

^ University of California Observatories/Lick Observatory, Uni- 
versity of California, Santa Cruz, CA 95064, USA 

George P. and Cynthia Woods Mitchell Institute for Funda- 
mental Physics and Astronomy, and Department of Physics and 
Astronomy, Texas A&M University, College Station, TX 77843- 
4242, USA 

Department of Astronomy, University of Massachusetts, 
Amherst, MA 01003, USA 

Department of Physics, Purdue University, 525 Northwestern 
Ave., West Lafayette, IN 47907, USA 

Department of Physics & Astronomy, University of California, 
Riverside, CA 92521, USA 

^'^ Deptment of Physics & Astronomy, University of California, 
Irvine, CA 92697 

Observatories of the Carnegie Institution for Science, 813 
Santa Barbara Street, Pasadena, CA 91101, USA 

Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, MS65, Cambridge, MA02138, USA 

^"^ Steward Observatory, University of Arizona, 933 N. Cherry 
St., Tucson, AZ 85721, USA 

^® Racah Institute of Physics, The Hebrew University of 
Jerusalem, Jerusalem 91904, Israel 



1. INTRODUCTION 

In the recent years, deep near-infrared imaging surveys 
have begun to yield a significant number of candidate 
galaxies at very high redshifts. The wide-field surveys 
from the ground have produced a handful of bright candi- 
dates at z w 7 (e.g., Ouchi et al. 2009; Hickey et al. 2010; 
Castellano et al. 2010; Capak et al. 2011; Hsieh et al. 
2012; Hathi et al. 2012), while the pencil-beam survey 
by the Hubble Space Telescope (HST) Wide Field Cam- 
era 3 (WFC3) within the historical Advanced Camera 
for Surveys (ACS) Hubble Ultra Deep Field (hereafter 
ACS-HUDF; Beckwith et al. 2006) and its two parallel 
fields (HUDF09; PI. Illingworth) have probed the fainter 
populations at z w 7-8 (Oesch et al. 2010; Bouwens et 
al. 2010, BIO; Bunker et al. 2010; McLure et al. 2010, 
MIO; Yan et al. 2010, YIO; Finkelstein et al. 2010) 
and possibly out to z « 10 (Yan et al. 2010; Wyithe et 
al. 2011; Bouwens et al. 2011a). The WFC3 Early Re- 
lease Science program (ERS, PI. O'Connell; Windhorst 
et al. 2011) has played an important role in connecting 
the "wide" and the "deep" ends of the exploration, al- 
lowing construction of z ss 7-8 samples at intermediate 
brightness levels (Wilkins et al. 2011; Bll; Lorenzoni 
et al. 2011, Lll). Meanwhile, the Hubble Infrared Pure 
Parallel Imaging Extragalactic Survey (HIPPIES; Yan 
et al. 2011) and the Brightest of Reionizing Galaxies 
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Survey (Trcnti ct al. 2011a) have been exploring the 
bright end of the population through a large number of 
random, discrete WFC3 pointings obtained during the 
HST parallel orbits. Surveys in foreground cluster fields 
utilizing gravitational lensing magnification have also re- 
sulted in a handful of promising candidates at 2 > 7-9 
(e.g., Richard ct al. 2006; Bradley et al. 2008, 2012; La- 
porte et al. 2011; Hall ct al. 2012), which complement 
the surveys in blank-sky fields. 

The Cosmic Assembly Near-Infrared Deep Extragalac- 
tic Legacy Survey (CANDELS; Pis: Faber & Ferguson; 
sec Grogin ct al. 2011 and Koekcmocr et al. 2011) 
employs a two-tiered, "wide and deep" strategy, which 
makes it uniquely positioned in bridging the significant 
gaps among existing surveys. In particular, its "Deep" 
component will cover ~ 125 arcmin^ in two fields upon 
completion, and its data will be ideally suited for study- 
ing galaxies at very high redshifts. In this paper, we 
report our preliminary results from the Deep observa- 
tions in the Great Observatories Origins Deep Survey 
(GOODS; Giavalisco et al. 2004) southern field, where 
wc use the nearly complete data set to study the galaxy 
population at z « 8. Wc use the following cosmologi- 
cal parameters throughout: — 0.27, I^a = 0.73 and 
Hq = 71 km s~^ Mpc~^. The quoted magnitudes are all 
in the AB system. 

2. DATA AND PHOTOMETRY 

As detailed in Grogin et al. (2011), the CAN- 
DELS/Deep region in the GOODS-S field is located in its 
middle stripe where we have the deepest GOODS Spitzer 
coverage, which also includes the ACS-HUDF. The data 
used in this study include 100% of the observations in 
F105W (hereafter Y105) and 83% of the observations in 
F125W ( J125) and F160W (iJieo) scheduled for this field. 
The typical exposure times per sky position are 8090, 
7450, and 7770 seconds in Y105, ^125 and iJieoj respec- 
tively. Our present analysis is confined to the region 
where the effective integration times are at least 4000 sec- 
onds in all three WFC3 bands, which is 62.9 arcmin^ in 
size. For the reason described below, we further exclude 
the region covered by the ACS-HUDF data, and the final 
effective area in our current study is 50.1 arcmin^. 

The products of the CANDELS data reduction (Koeke- 
moer et al. 2011) include the science mosaics and the as- 
sociated "root mean square" (RMS) maps. These RMS 
maps account for correlations between pixels that are in- 
troduced during the data reduction procedures (see e.g., 
Dickinson et al. 2004). Such correlations tend to reduce 
the measured pixel-to-pixel background noise in the HST 
images, and can lead photometry programs to underes- 
timate the actual photometric uncertainties. We avoid 
this problem by computing uncertainties using the nor- 
malized RMS maps, summing their values in quadrature 
over the measurement apertures, and including an addi- 
tional term for Poisson noise from the source fiux (gen- 
erally negligible for very faint dropout galaxies). The 
values of the CANDELS RMS maps have been validated 
both by comparison with autocorrelation measurements 
of the image noise, and by analysis of the dispersion of 
measurements in apertures randomly distributed within 
empty regions of the images. Our final mosaics have a 
pixel scale of 0'.'06, and are tied to the same astrometric 
frame as the ACS data in GOODS. 



The GOODS ACS data are used in our study to con- 
firm the optical non-detections of the candidates (see 
§3.1). These data have been reprocessed to make mosaics 
at the same 0f.'06 pixel scale for CANDELS. On average, 
the 5-a sensitivities (as measured from the RMS maps) 
within a O'.'2-radius aperture in ACS F435W, F606W, 
F775W, F850LP (hereafter B435, Veoe, »775, zs^o) and 
WFC3 Yios, J125, and Hieo are 28.04, 28.25, 27.65, 27.48, 
28.15, 28.05, and 27.82 mag, respectively. 

The HUDF region presents a complexity in our analy- 
sis. Both the HUDF09 WFC3 IR data and the ACS- 
HUDF data are much deeper than their counterparts 
that cover the rest (and much wider) area of interest. 
If we were to use these deeper data in our candidate se- 
lection, it would be very difficult to uniformly carry out 
the simulations that are necessary for the follow-up anal- 
ysis (sec §3.4). For simplicity, wc therefore exclude the 
ACS-HUDF region from this current work. 

In addition to the GOODS ACS data, we also have 
new ACS data in F814W (hereafter hu) that CANDELS 
took in coordinated parallel mode to the WFC3 IR ob- 
servations (sec Grogin et al. 2011). These new ACS 
data suffer severely from the degraded charge transfer 
efficiency (CTE) of the aged CCDs, however, and the 
empirical correction of Anderson & Bedin (2010) were 
applied first before they were processed by the CAN- 
DELS reduction pipeline. 

The photometry in WFC3 and ACS was done together 
by running the SExtractor software (Bertin & Arnouts 
1996) in dual-image mode. We used the sum of the 
J125 and the i?i6o mosaics, each weighted by the inverse 
square of their RMS maps, as the detection image. A 
3-pixel FWHM, 7x7 Gaussian filter was used to con- 
volve the detection image, the threshold was set to 0.8- 
fT, and at least 4 connected pixels above this threshold 
are required. SExtractor provides a number of meth- 
ods to measure flux, among which MAG_APER, MAG_ISO 
and MAG_AUTO are commonly used in dropout selections. 
MAG_APER magnitudes are calculated within a circular 
aperture of fixed size. As most high-redshift galaxies 
are compact, using these magnitudes can still be appro- 
priate when the circular aperture size is carefully chosen 
so that it maximizes the S/N and minimizes the loss of 
light for the majority of the sources of interest. While 
it has the advantage of being robust in extracting very 
faint sources, MAG_APER has the disadvantage that real 
galaxies, which have different light profiles, have differ- 
ent fractions of light lost for a fixed apertm-c. MAG_ISO 
magnitudes arc calculated within isophotal areas, which 
are set by the detection threshold and hence trace the dif- 
ferent shapes of real galaxies. For this reason, MAG_ISO 
usually can be tuned to maximize S/N for galaxies of 
different morphologies. However, it docs not capture 
the total light of a galaxy. MAG_AUTQ magnitudes arc 
calculated using adaptive elliptical apertures based on 
the first moment of the light distributions of galaxies 
(Kron 1980), and are the most suitable for deriving total 
magnitudes if an appropriate "Kron factor" is chosen. 
However, using a large Kron factor to capture total light 
usually means that the apertures are set too large for 
very faint sources such that their S/N are not optimized. 
Bearing both the advantages and the caveats in mind, 
we adopted the MAG_AUTO magnitudes (using the default 
[Kron factor, minimum radius] = [2.5, 3.5]) in this work. 
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Dropout selection uses hard boundaries in color space 
(see §3.1). As a result, photometric errors cause objects 
be scattered in and out of the selection area: an ob- 
ject whose intrinsic colors meet the selection threshold 
could be missed because of the errors in its measured 
colors, and vice versa. While this effect can be taken 
into account statistically by applying the correction for 
the "effective volume" (see §3.4), the candidate samples 
constructed using different photometry could differ sig- 
nificantly on an object-by-object basis. To demonstrate 
this effect, we also constructed a separate sample using 
MAG_ISO magnitudes, which is presented in the appendix 
and is compared to the MAG_AUTO sample. In both cases, 
we only kept the sources that have S/N > 5.0 in at least 
one band among J125 and ilieO) measured within the 
MAG_AUTO or the MAG_ISO apertures. 

We also make extensive use of the GOODS Spitzer 
IRAC data (Dickinson et al. 2004), which also have 
the same sky projection as other data here, and have a 
pixel scale of C.'S. We mostly concentrate on the 3.6 and 
4.5 /xm channels (hereafter [3.6] and [4.5]) in this study. 
We used a version of the GOODS photometry where the 
detection was done with a "mcxhat" filter to optimize 
de-blending. We adopted MAG_APER within a l'.'5-radius 
(i.e., 2.5 pixels) aperture, and obtained the "total mag- 
nitudes" by applying the corrections as determined in 
the GOODS program where the aperture of this size has 
been used (e.g., Yan et al. 2005), which are —0.55 and 
—0.60 mag for [3.6] and [4.5], respectively. 

3. yio5-DROPOUT SELECTION OF CANDIDATE 
GALAXIES AT Z w 8 

In this section, we describe our samples of candidate 
galaxies at z « 8 selected as Yios-dropouts and the pos- 
sible biases in our selections. 

3.1. Color Criteria 

Our first color criterion for Fio5-dropout selection is 
^105 — J125 > 0.80 mag. This is the same as that em- 
ployed by BIO and YIO, and is significantly larger than 
that of Y105 — J125 > 0.45 mag used in Bll and Oesch 
et al. (2012) where those authors have opted to included 
more galaxies at lower redshifts (starting at 2: > 7.2) in 
their Yios-dropout selection. Our larger color decrement 
threshold is sensitive to Lyman-break at 2; > 7.7. When 
the objects have S/N < 1 in Y105, we use their 2 cr up- 
per limits and only keep those that have Y105 — J125 > 
0.80 mag as calculated using these conservative; upper 
limits. Our second criterion is that a Yios-dropout should 
have S/N < 2 in all four bands of GOODS ACS images 
{Bi35VeQein5Zs5o). This is to implement the invisibility 
requirement in these "veto" images, which is a neces- 
sary condition (though not sufficient) that the detected 
^105— "^125 color decrements in our candidates are not fea- 
tures in galaxies at low redshifts (for example, the 4000 A 
break at z « 2) that could mimic the Lyman-break, be- 
cause such low-redshift galaxies could have detectable 
fluxes in these bands. Setting the threshold to a quanti- 
tative value of S/N < 2 is somewhat arbitrary, however. 
Based on experience, adopting S/N < 2 for invisibility 
usually is a good compromise between minimizing the 
incompleteness and keeping the size of the initial candi- 
date sample manageable for the later visual inspection 



(see below). To exclude as much as possible the contam- 
ination from the red galaxy population at low redshifts, 
we impose a third criterion of J125 — -ffieo < 0.3 mag. 
This also largely limits our selection window to 2 < 8.7, 
and could be biased against genuine galaxies at z « 8 
that have unusually large dust reddening or extraordi- 
narily old stellar populations should they exist. Figure 1 
shows our selection in the yios ^125-^160 space. 

Each of these criteria introduces incompleteness, for 
which we will correct statistically (see §3.4). 

3.2. Visual Inspection 

The initial candidates have been visually inspected to 
ensure that they are legitimate sources in the IR images 
and that there is no reason to suspect that their pho- 
tometry might be unreliable. Examples of problems that 
are not accounted for in the statistical noise model (and 
thus require visual inspection) include diffraction spikes, 
poorly rejected cosmic rays, inappropriate deblending, 
and bad background subtraction due to nearby neigh- 
bors. While such issues are rare in the entire photomet- 
ric catalog, they can be selected when looking for objects 
with unusual colors, such as dropouts. Visual inspection 
step is necessary to reduce the rate of contamination, 
and is commonly adopted in dropout selections by vari- 
ous groups. It is somewhat subjective; different inspec- 
tors might not agree on a particular candidate, and even 
the same inspector might not be able to completely re- 
produce at different time his/her results on a large set of 
samples. Nevertheless, visual inspection is irreplaceable 
to get a reliable sample. 

Ideally, the simulation discussed in §3.4 should also 
include this step when deriving the incompleteness cor- 
rection. However, it is impractical to examine the 90,000 
artificial galaxies inserted during the simulation. Never- 
theless, we do not think that many of these would have 
been inappropriately rejected by visual inspectors, be- 
cause the inserted galaxies do not resemble the artifacts 
or those with skewed photometry that we typically re- 
ject. Relative to the 90,000 inserted objects, we suspect 
only a handful would have these problems. Thus we do 
not think that the failure to visually inspect the artifi- 
cial galaxies has a significant effect on our computation 
of the incompleteness and the effective volume. 

In our visual inspection process, the initial candidates 
were first examined by one of us (HY) to get rid of the 
most obvious contaminators (> 300 objects), such as 
image defects close to the field edges, diffraction spikes 
of bright stars, bad pixels, residuals of IR persistence, 
etc. The remaining candidates (< 200) were passed to 
six persons for independent inspections using a three- 
grade system of "firm acceptance" , "firm rejection" and 
"border-line case" . To be included in the final samples, 
a candidate should be firmly accepted by at least three 
inspectors, and should not be firmly rejected by more 
than one inspector. As described below, this approach 
has allowed us to study the possible systematic biases 
caused by visual inspections. 

We note that the contemporaneous CANDELS /814 im- 
ages have also been used in our current selection process. 
However, the flux upper limit estimate is subtle in these 
CTE-restored images when a source is undetected in sin- 
gle exposures, and a rigorous treatment to be adopted 
by our team is still yet to be finalized. To avoid possi- 
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ble changes to our sample in the future because of this 
reason, at this stage we refrain from directly using the 
/814 upper limits or S/N values in the quantitative veto 
process, and just use the /8i4 images in the visual inspec- 
tion step (nevertheless, we do confirm that our current 
final candidates indeed formally satisfy the invisibility 
requirement in /814; see §3.3). 

In addition to the usual examination of the mosaics, 
we have taken extra steps for assurance. In the IR, wc 
have followed Yan et al. (2011) and examined the fi- 
nal candidates in their single-epoch WFC3 science im- 
ages and the associated "data-quality" flags to confirm 
that they are not due to instrument defects. In the opti- 
cal, we have constructed the "x^" sum image of the four 
GOODS ACS B435V606«775Z850 bands for our candidates. 
Based on Szalay et al. (1999), the sum of these im- 
ages optimally combines the signals from different bands 
and extend to a fainter limit that achieved by individual 
bands. The x^ images of our final candidates all have 
X^ < 3.73, which is consistent with being drawn from a 
random sky distribution (Szalay et al. 1999) and thus 
further confirms their non-detections in the ACS. 

3.3. Discussion of Sample 

The final candidates in our sample are listed in Table 
1. Their images are displayed in Figure 2. 

3.3.1. "Negative Image" Test 

To investigate the possible contamination of our sam- 
ple by noise spikes due to background fluctuations, wc 
carried out the so-called negative image test (e.g., Dick- 
inson et al. 2004). The image mosaics described in §2, 
which all have zero mean background by construction 
(see Koekemoer et al. 2011), were multiplied by —1 to 
make any positive pixels to negative and vise versa. SEx- 
tractor was run in dual-image mode again on these nega- 
tive images in the same way as in §2, using the weighted 
sum of the the negative J125 and Hica images as the 
detection image. After applying the same selection cri- 
teria for our lio5-dropouts as in §3.1, only three "ob- 
jects" survived, all of which have J125 > 27.5 mag, and 
all only show up in the negative J125 image. Their mor- 
phologies have odd shapes and are typical of noise spikes, 
and would be easily rejected during the visual inspection 
step. Therefore, we conclude that the contamination due 
to noise spikes is negligible in this work. 

3.3.2. Interlopers 

Our stringent requirements of non-detection in the 
ACS (individual bands and x^ sum) as well as J125 — 
^^160 l£ 0.3 mag arc effective in removing contaminations 
from red galaxies at intermediate redshifts (sec Figure 1). 
In addition, the IRAC data for the isolated candidates 
support that they are unlikely contaminators of this kind, 
because such objects, if very weak in optical, should have 
very red colors at 1-4 fim such that J125 — [3.6] > 2 mag 
(e.g., Yan ct al. 2004). Our candidates that have IRAC 
flux measurements are all at J125 — [3.6] ^1.3 mag. 

Another possible contaminant could be cool Galactic 
dwarf stars. Our criterion of Yios — J125 > 0.80 mag is 
large enough to exclude most such contaminations. In 
addition, the number density of M/L/T-dwarfs in a high 
galactic latitude region such as the GOODS-S is likely 



negligible (Ryan et al. 2011). Furthermore, our detailed 
morphological analysis in the accompanying paper (R. 
Ryan et al. 2012, ApJ submitted) shows that our bright- 
est candidates are more consistent with extended light 
profiles than point-like. 

Using the 4 Ms Chandra catalog of Xue et al. (2011), 
we find no X-ray counterparts to our candidates within 
6 . There is thus no evidence that these objects are AGN 
rather than normal galaxies. 

3.3.3. Comparison to the Sample from Oesch et al. (2012) 

Another paper based on the similar data set as ours 
was submitted by an independent group right after ours 
(Oesch et al. 2012, hereafter 012), and here we compare 
our sample to theirs based on the latest update available 
to us (P. Oesch, private communication, 2012). 

In the overlapping region, the WFC3 IR data are essen- 
tially the same in both studies. However, 012's reduction 
of the CANDELS WFC3 data is different from ours, and 
hence the end products of the science mosaics and the 
associated RMS maps are different. The photometry is 
also different. 012's detection image is a x^ image de- 
rived from the J125 and i?i60 images, whereas ours is the 
weighted sum of the J125 and iJieo- While they used 
MAG_AUTO as we also did, they adopted a smaller aper- 
ture (1.2 Kron) to compute colors. In addition, they 
stacked all the available ACS data in this region. We 
refrained from combining the late ACS data that suffer 
from the CTE problem with the earher GOODS V2.0 
-8435^06^775-2850 ACS data, and for the same reason we 
only used the /814 data in the visual inspection step (see 
§3.2). Finally, 012 adopted rather generous color crite- 
ria of yio5 — J125 > 0.45 mag and J125 — -ffieo < 0.5 mag, 
while we used Y105 — J125 > 0.80 mag and J125 — i?i60 < 
0.3 mag. 

All these factors contribute to the differences between 
the 012 sample and ours. It is interesting and instruc- 
tive to compare the inclusion of the 012 candidates in 
our samples, and vice versa. Based on the photome- 
try reported in 012, nine of their 14 yio5-dropouts in 
the same field would be expected to enter our samples. 
Our sample, as shown in Table 1, include five of them. 
Among the four candidates not in our sample, three of 
them, CANDY-2209246371, 2432246169 and 2277945141 
have S/N < 5 in J125 and H\qq and are not included 
in our MAGJVUTO catalog, and the other one, CANDY- 
2272447364, has Yios - J125 = 0.18 based on our pho- 
tometry. 

We note that one of our brightest candidates, 
AUT0_100, is not in the 012 sample. 012 beheves that 
this object has a 2.4 a detection in /8i4, derived in the 
conventional way using SExtractor. While our selec- 
tion does not use the formal S/N estimates in /814 (see 
§3.2), we have also derived such estimates and confirmed 
after-the-fact that all our final candidates formally have 
S/N< 2. This particular object has S/N< 2 in both the 
MAG_AUTO and the MAG_ISO apertures (see also appendix) . 
We attribute this discrepancy to the differences in pho- 
tometry as outlined in §3.3.2. 

To conclude, the level of overlap and discrepancy be- 
tween these two samples is all well expected given the 
level of differences in their constructions. We further 
demonstrate this in the appendix. 
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3.4. Redshift Selection Functions 

A given choice of Lyman-break color criteria, applied 
to a particular data set, leads to a probability distribu- 
tion as a function of redshift and magnitude, P(m, z) 
for selecting galaxies as dropout candidates. This prob- 
ability distribution, in turn, corresponds to an effective 
volume (Veff) for Lyman-break selection as a function of 
magnitude (Steidel et al. 1999). We have carried out ex- 
tensive simulations to derive P(m, z) and Veff, in order 
to correct for the sample incompleteness that is inherent 
in the dropout technique. The results are the redshift 
selection functions at different brightness levels, which 
are shown in Figure 3. 

To derive these functions, we put 90,000 artificial 
galaxies into the WFC3 and the ACS mosaics and calcu- 
lated their rate of recovery P{m, z) after applying the 
same color criteria as we did when constructing our 
dropout samples. These artificial galaxies have differ- 
ent fluxes and morphologies, with 70% being exponential 
disks and 30% having the de Vaucouleurs profile. Their 
sizes follow a log-normal distribution and 9{f) oc f'^'^^, 
where 9 is the peak size at flux / and is normalized to 
0'.'2 at HiQQ = 26 mag. The input spectrum is from 
the models of Bruzual & Chariot (2003; BC03), and 
has a constant star formation history (SFH) and an age 
of 50 Myr. We adopted the extinction law of Calzetti 
(2001), and assumed that the distribution of E(B-V) is 
a Gaussian function with mean (E(B-V)) — 0.1 mag and 
scatter a = 0.1 mag, restricted to E(B-V) > 0. The SEDs 
of the simulated galaxies were attenuated by the line-of- 
slight H I absorption according to the recipe of Meiksin 
(2006), which is very close to that of Madau (1995). 

By calculating the ratio of the number of the recovered 
galaxies and that of the input ones, we obtained the selec- 
tion function for the incompleteness correction. Figure 
3 presents it as a function of J125 MAG_AUTD magnitudes, 
as these magnitudes are taken as the total magnitudes, 
and are the closest to the input total magnitudes of the 
simulated galaxies. 

3.5. IRAC Counterparts 

We have searched for IRAC counterparts of our Y105- 
dropouts within r = 0'.'6, and the results are also sum- 
marized in Tabic 1. Only two objects are securely de- 
tected in IRAC, namely, AUT0_035 and AUT0_293. As 
we will discuss later, their WFC3 versus IRAC colors 
are very consistent with being at high redshift. Objects 
AUT0_204, 212, 368 are in isolated regions, and they 
are invisible in the GOODS IRAC data. For all these 
objects, their 2-a flux upper limits in [3.6] and [4.5] chan- 
nels are calculated within a r = 1'.'5 aperture. 

All other objects in our samples are either blended 
with (denoted by "B" in Table 1) or severely contam- 
inated ( "C" ) by foreground neighbors in IRAC and thus 
no useful information can be obtained. 

3.6. Summary of Sample Limitations 

Our sample suffer from the same kinds of incomplete- 
ness and contamination that are intrinsic to any dropout 
samples because of their selection technique. While a 
number of these issues have been mentioned in the above 
sections, we review them here before proceeding to dis- 
cuss the statistical implications of our samples to astro- 
physical questions. 




-1 1 

"'las ~ ^160 

Fig. 1. — Selection of Yios-dropouts in the Yi05^l25J^i60 color- 
color diagram. The small dots are field objects that have S/N< 2 
in ACS. The grey area is our selection region, where the brighter 
candidates (J125 < 26.2 mag) are plotted as red filled squares, and 
the fainter candidates as red open squares. The blue-magenta track 
shows a young galaxy (using a BC03 model of age 100 Myr and no 
dust extinction) a.t z = 5—10, with the blue section for z < 7.7 and 
the magenta section for z > 7.7. The orange stars show the colors 
of Galactic brown dwarfs (Leggett et al. 2002). The green symbols 
show the colors of a typical red galaxy at 2 1-3 simulated using 
a BC03 model (r = 50 Myr and age of 2.0 Gyr). 
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Fig. 2. — Images of the Yios-dropouts in the AUTO sample, 
each S'.'l on a side. The locations of the candidates are at the 
center marked by red circles, which are 0'.'5 in radius. The com- 
mon objects in the comparison ISO sample (see appendix) are also 
labeled. 



There are two main factors that contribute to the in- 
completeness, which cause significant differences among 
samples as demonstrated in §3.3.1 and 3.3.2. First of 
all, an object whose intrinsic brightness is above our sen- 
sitivity limit could be undetected because its observed 
brightness could fall below the detection threshold due 
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Fig. 3. — Redshift selection functions derived from simulations 
of our color selection in different magnitude ranges. The labeled 
magnitude ranges are based on MAGJHJTQ in J125. 



to noise. Secondly, a genuine 2; « 8 galaxy that is de- 
tected in source extraction could still escape our color 
selection because it could be "scattered" out of the color 
selection region (see Figure 1) due to photometric error, 
or it could be incorrectly "vetoed" because of some unre- 
lated events within its vicinity (for example, noise spikes) 
that accidentally increase its measured S/N in the veto 
images above the adopted quantitative veto threshold. 
Both effects have been taken into account by our simu- 
lations, and therefore can be corrected statistically. 

The major factor that causes the contamination is the 
inclusion of low-redshift objects with spectral features 
that could mimic the Lyman-break at z « 8. Such con- 
taminators could be selected because their Y105 — J125 
colors meet our color threshold, either intrinsically, or af- 
ter being boosted by photometric error (i.e., being scat- 
tered into the selection region). Using the i/ieo data, 
our secondary J125 — iJieo criterion reduces the number 
of contaminators by rejecting those that have red SEDs. 
For other contaminators, the best way to reject them is 
to check their visibility in the veto images. As discussed 
earlier, the quantitative criterion of S/N < 2 gets rid of 
most of such objects, and the final step is the visual in- 
spection. The caveat here is that we implicitly assume 
that the veto images are deep enough in detecting the 
contaminators in the blue bands, but this is not guaran- 
teed. Currently we do not have a reliable way to assess 
this effect. 

We point out that the rejection of contaminators based 
on their visibility in the veto images also introduces in- 
completeness, and that this happens both in the quanti- 
tative pre-screening (requiring S/N < 2) step and in the 
last visual inspection step. While the former can be sta- 
tistically corrected by the simulations, the latter cannot. 

4. DISCUSSION 



Our samples are unique in two aspects. (1) We have 
two very bright candidates at J125 < 26.0 mag, which 
are at least ~ 1 mag brighter than the brightest Y105- 
dropouts previously found. They are comparable in 
brightness to the brightest Yogs-dropouts reported by 
Yan et al. (2011) and Trenti et al. (2011a,b) in their 
WFC3 parallel surveys, but could be at higher red- 
shifts because the redshift window of this current Y105- 
dropout selection is higher than that of the existing yogs- 
dropouts. (2) Two other Yio5-dropouts, which are not 
among the brightest in the WFC3 bands, are securely 
detected in IRAC, which is in sharp contrast to previ- 
ous results where only non-detection in IRAC have been 
reported for Yios-dropouts. Here we discuss the implica- 
tions of both. 

4.1. Stellar Population and Stellar Mass 

We investigate the stellar populations of the IRAC- 
detected Fios-dropouts, AUT0_035 and AUT0_293, by ana- 
lyzing their spectral energy distributions (SED) through 
template fitting. The approach here is similar to that 
of Finkelstein et al. (2011). We first use the EAZY 
code (Brammer et al. 2008) and the photometry of 
the GOODS ACS & CANDELS WFC3 data to estimate 
photometric redshifts {zphot)- The full SED (including 
the IRAC photometry) is then fitted to a suite of tem- 
plates based on the updated models of BC03 (the so- 
called "CB07" models). A Salpeter initial mass func- 
tion (Salpeter 1955) and the metallicities of 0.02-1^0 
are adopted. The templates have a range of exponen- 
tially decreasing and increasing SFHs (see also Papovich 
et al. 2011), and can include nebular lines based on the 
number of ionizing photons and metallicity of a given 
model (B. Salmon et al. in preparation). We assume 
the extinction law of Calzetti (2001) with E{B - V)=Q- 
0.5 mag and the H I absorption as formulated in Madau 
(1995). The results are summarized in Figure 4. We 
note that fitting redshift and other properties simultane- 
ously does not signficantly change these results, and the 
differences are captured in the errors that we quote here. 

The best-fit photometric redshifts are 8.5I45 

8.910.5 fo'" AUT0_035 and 293, respectively, which are 
consistent with the redshift window of our color selec- 
tion. Adding the formal /814 limit to the fitting pro- 
cess does not change the results. We fit the SEDs 
with and without the contribution from nebular lines, 
and both results suggest high stellar mass (A^*) in 
the range of 10^ to 10^°Mq for both objects. Includ- 
ing nebular emission lines, we obtain 2.btil x and 
0.9j:o ° X 10^ Mq for AUT0_035 and 293 respectively, while 
without nebular emission lines these are 10.5j:8:^ X 10^ 
and 2.6lJ'4 ^ -^0^ respectively. The error bars reflect 
the 68% confidence level of the fit. Labbe et al. (2010a) 
derived the stellar masses of galaxies at lower redshifts 
oi z K, 1 using models without nebular emission, and 
the average value is ~ 10^ Mq. The two objects in our 
sample have comparable or even larger values using sim- 
ilar models, which are surprisingly high at such an early 
epoch. 

Using the M.^ estimates, we obtain the first measure- 
ment of the mass function (MF) of galaxies at z 8 
at the high-mass end. As we only have two objects 
whose values are reasonably similar, we opt to 
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count them within one mass bin of a ±0.4 dex bin size. 
The two different sets of models result in significantly 
different mass estimates, and therefore the bin center 
is different for each case. Assuming a top-hat selec- 
tion function of our survey within 7.7 < z < 8.7, we 
get Hlog§^)\9.2±oA = (1.9l^:^) X 10"' Mpc'^ dex'i 
if we adopt the values derived using the models with 
nebular emission. If we apply the effective volume 
(Veff) correction of the survey instead (see §3.4), we 
obtain (S.lti!^) x 10"^ Mpc'^ dex'^. Here we use 
^eff = fVeff{mj)dm,j = J f dmjdzP{m,j,z)dV/dz, 
where dzdV/dz is the unit co- moving volume at red- 
shift z, and P{mj,z) is the redshift selection function 
at different magnitudes mj as derived through simula- 
tions (see Figure 3). The error bars here reflect the 68% 
interval of the uncertainties caused by the Poisson noise 
in the sample. For the case of using the TW* estimates 
based on the models without nebular lines, all the above 
values are applicable at </'(Zo(7-^)|9.7±o.4- 

The contribution of these objects to the global stellar 
mass density, p,, at z w 8 within our survey volume is 
i^.Stli) X 10^ and (12.61:^;^) x 10^ Mq Mpc'^ based 
on the models with and without nebular emission lines, 
respectively, if assuming a top-hat selection function at 
7.7 < z < 8.7. The errors reflect the 68% confldence level 
of the fitting results. If we apply the same correction 
for Veff as above, we obtain — (5.5^2^4^) ^ ^^'^ 
(21.21^5^3) X 10-* Mq Mpc-3. 

We derive the mass-to-light ratio of these two objects 
using their 4.5 /zm flux as a proxy to Ly- Using the 
fit results with contribution from nebular lines, we ob- 
tain M*/Lv = 0.071^5 ^4 and OMtom for AUTD_035 and 
293, respectively. Without the contribution from nebu- 
lar lines, we obtain M^/Ly = 0.3lt[;;^^ and 0.12l[J;°^. 
The latter values are in general agreement with those ob- 
tained by Labbe et al. (2010a, b) through stacking anal- 
ysis of the IRAC data of the Yogs-dropouts and the l^os- 
dropouts in the WFC3 ERS field (see also Gonzalez et al. 
2011) and the HUDF using similar models without neb- 
ular lines. On the other hand, Labbe et al. (2010b) esti- 
mated that nebular lines have a small impact on A^*/Ly, 
reducing it by ~ 0.2 dex. Our results suggest a stronger 
effect, reducing A4*/Lv by ^ 0.5-0.6 dex. 

4.2. Constraint on the LF 

Here we use our sample to constrain the LF at z w 8. 
We adopt the J125 MAG_AUTO values of the candidates as 
their total magnitudes. 

After applying the correction for Veff, we 
obtain the stepwise LF in the five 0.5-mag 
bins (25.45, 25.95, 26.45, 26.95, 27.45) mag as 
= (2.5t^o:«,2.5t^J,2.8t«:i3.6t«i,24.2™) x 
10 ^ Mpc ^ mag ^. These number densities are shown 
in Figure 5, and compared to the predictions from a 
number of Schechter LF estimates at z ~ 8 (MIO; YIO; 
Bll; Lll). The black squares are our observed densities, 
while the red squares are the densities after correcting 
for Veff. The error bars represent a Bayesian 68% 
credible interval, indicating the central 68 percentile 
range for the posterior distribution of the true number 
density assuming Poisson statistics. These uncertainties 
account for the incompleteness, but do not account for 
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Fig. 4. — Summary of SED fitting results for AUT0_035 (left) and 
293 (right). The results obtained from the fit with and without 
the contributions from nebular lines are coded in blue and red, 
respectively. The top panels show the observed SEDs (grey circles) 
and the best-fit templates, and the insets display their 3.6 and 
4.5 lira image (18'.'6x 18'.'6). The bottom panels show the likelihood 
functions of Zphot S'Ud , and the best-fit values are indicated by 
the vertical lines. 

any systematic uncertainties from the contamination 
due to low-z interlopers. For ease of direct comparison 
to observations, surface density and apparent magnitude 
scales are also provided on the same figures to present 
these results in terms of differential number densities 
versus apparent J125 magnitudes. For comparision, we 
also plot the observed densities (i.e., before corrections 
for their corresponding Veff values) extracted from the 
sample of YIO in the HUDF09 proper and those of 
Bll in the HUDF09, HUDF09P1 and HUDF09P2 after 
applying the additional criteria of Y105 — J125 > 0.8 mag 
and J125 - Hieo < 0.3 mag. 

An intriguing feature that our samples reveal is that 
there could be an excess of bright objects at J125 < 
26.2 mag with respect to any of the existing Schechter LF 
estimates at z « 8 in the literature would predict. The 
inferred number density at this bright end is even higher 
than the prediction from a non-evolution z ss 7 LF (e.g., 
the one from Bll). This excess is still evident with the 
observed counts before applying the incompleteness cor- 
rection. The observed counts at J125 < 26.2 mag in Fig- 
ure 5 comes from objects AUTD_048 and AUT0_100. Based 
on their photometry in J125 and assuming z = 8, they 
have Muv of —21.54 and —21.21, respectively. While 
such high luminosities have been observed in the spectro- 
scopically confirmed z w 7 galaxies of Ono et al. (2012), 
they are 5-9 x more luminous than any previously known 
Yi05-dropouts. Such an excess is unlikely due to cosmic 
variance, whose effect is expected to be small compared 
to the Poisson errors for our observed number density if 
we use the formalism of Trenti & Stiavelli (2008). 

We note that, however, that 012 do not see such an 
excess, even though their sample also includes our bright- 
est candidate AUTQ_048. The main reason for this dis- 
crepancy could be that 012 use a larger volume in the 
calculation. They have included part of the CANDELS 
Wide field in the GOODS-S region, which is a major 
contributor to their larger volume. However, this area 
does not have sufficiently deep Y105 data to apply our 
selection criterion of F105 — 0^125 > 0.80 mag down to 
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Fig. 5. — Constraints on the very bright-end of the LF at 2 ^ 8 
based on the sample of Yios-dropouts in CANDELS/Deep from 
this current work. The black squares are the differential surface 
densities inferred from the number counts in our sample, while 
the red squares are those after correction for V^ff with respect 
to the volume within 7.7 < z < 8.7. The count predictions over 
7.7 < z < 8.7 are from the Schechter LF estimates at 2 8 
from Bll, MIO, YIO, and Lll are shown as various black curves, 
together with a non-evolution one at z ^ 7 from Bll as the blue 
curve. For comparison, different symbols at the faint-end show 
the raw (observed) surface densities of lios -dropouts based on the 
following samples: YIO (HUDF - stars), Bll (HUDF - crosses; 
HUDFPl - pentagons; HUDFP2 - triangles). 

Ji25 ~ 26.0 mag. 012 have adopted a less stringent cri- 
terion of Yio5 — Ji25 > 0.45 mag, which allows them to 
construct a Yios-dropout sample in this area at the cost 
of including candidates at lower redshifts {z « 7.2). The 
advantage of their approach is that they can incorporate 
the WFC3 ERS field in the GOODS-S region as well, 
which further increases their survey volume. While the 
WFC3 ERS field has Fogs but not Y105 data, the redshift 
range probed by their Iggs-dropouts in this area overlaps 
that of their Yio5-dropouts elsewhere, and therefore 012 
has combined the Yogs" ^-nd the Yios-dropouts to address 
the LF averaged over a wide redshift range (z w 7-8.7). 

In contrast, our color criteria select galaxies in a higher 
range of redshifts, 7.7 ^ z < 8.7, and the bright-end 
excess (with respect to a smooth Schechter function) that 
we find suggests that our understanding of the LF at 
z w 8 is still far from being complete. The significance of 
this excess is still subject to small number statistics, and 
we cannot rule out the possibility that one or both of our 
brightest objects could be interlopers at low redshifts. 
However, we should point out that similar bright-end 
excess at z ~ 8 has also been suggested in other studies, 
for example in the bright Fogs-dropout study in HIPPIES 
(Yan et al. 2011) at similar depths, in the ground-based 
survey of Laporte et al. (2011) at a brighter level, and 
in the much wider survey in Capak et al. (2011) at an 
even brighter limit 

22 While Bowler et al. (2012) derived different Zph i~ 2) for the 



4.3. Implications for Star Formation in the Early 
Universe 

Our z « 8 candidate galaxy samples have two very 
luminous objects and two high stellar mass objects. We 
have presented tentative evidence that the bright end of 
the galaxy LF at z ss 8 might not follow the exponential 
cut-off of the Schechter function. While the stellar mass 
function at this redshift has not yet been determined, we 
have provided the first, albeit still tentative measurement 
at the high-mass end and shown that some galaxies with 
stellar masses as high as a few x 10^ M© might already 
in place at z w 8. 

The physical origin of the exponential cut-off at the 
bright end of the LF or the high mass end of the MF 
seen at lower redshifts remains a matter of debate. At 
z ~ 0, star formation is "quenched" in halos with masses 
greater than ~ 10^^ Mq ("quenching mass"), perhaps by 
AGN feedback (e.g., Croton et al. 2006; Somerville et al. 
2008; Gabor et al. 2011). However, at higher redshift 
(z < 2), the quenching mass may be higher (Dekel et 
al. 2009; Behroozi et al. 2010), and the existence of 
massive, rapidly star forming galaxies at 2 > 2 is well 
established. At z ~ 2, the observed exponential cut- 
off in the rest-frame UV LF appears to be due to dust 
(Reddy et al. 2010) rather than quenching. At very high 
redshifts, it has been noted before that the very blue 
rest-frame UV colors of z w 7 candidate galaxies (some 
of them confirmed) suggest that they may contain little 
dust (e.g., Oesch et al. 2010; M10;'y10; Bunker et al. 
2010). Thus, at z « 8, when (1) the quenching mass 
is much higher than the typical halo mass, and (2) dust 
has little effect, perhaps we should not expect to see an 
exponential cutoff in the luminosity or mass functions of 
UV-bright galaxies. 

It is interesting that the estimated stellar masses and 
number densities could imply a rather high efficiency of 
conversion of baryons into stars. For the lower stellar 
mass estimates, which are obtained using the models 
with nebular emission (the average is ^ 2 x 10^ Mq), 
and assuming that ^ 20% of the available baryons have 
been converted to stars, the implied halo masses are on 
the order of ^ 10"'^^ M©, for which the expected num- 
ber density of dark matter halos in the currently favored 
LCDM cosmology is comparable to the observed num- 
ber density of objects. This is consistent with the pre- 
dicted star formation efficiencies and host halo masses 
at z w 8 from cosmological hydrodynamic simulations. 
However, if the true stellar masses are higher, or the star 
formation efficiency is lower, a rapidly growing tension 
arises between the number density of dark matter halos 
and the observed number density of galaxies (above halo 
mass ~ 10^^ Mq, the halo number density declines by 
about two orders of magnitude for a factor of 3 in- 
crease in mass). While converting 20% of the available 
baryons into stars may not sound excessive, this is in 
fact the maximum value that has been inferred at any 
epoch. Due to the presumably very low metallicity of 
the gas in these early objects, we might have expected 
much lower star formation efficiencies than are seen lo- 
cally (e.g. Krumholz & Dekel 2010). 

candidates in Capak et al. (2011), the true nature of these objects 
still remains uncertain. 
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5. SUMMARY 

In this work, we search for candidate galaxies at ^; « 8 

in the CANDELS Deep GOODS-S field and study their 
properties. Our sample of yio5-dropouts hints that the 
number density of z « 8 galaxies at the bright-end might 
be higher than expected from the previous Schechter LF 
estimates, which lends support to the suggestion made by 
a number of earlier studies that there could be a bright- 
end excess in the galaxy number density at very high red- 
shifts. Furthermore, two of our candidates are securely 
detected at 3.6 and 4.5 iira in Spitzer data. These are the 
first lios-dropouts individually detected at these wave- 
lengths. Their derived stellar masses are on the order of 
~ 10^ Mq, from which we obtain the first measurement 
of the high-mass end of the galaxy stellar mass function 
at 2; « 8. If the high number densities of very luminous 



and very massive galaxies at z « 8 are real, they could 
imply a large stcUar-to-halo mass ratio and an efficient 
conversion of baryons to stars at very early time in the 
cosmic stellar mass assembly history. 
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TABLE 1 

MAG_AUTO-SELECTED Vios-DROPOUTS IN CANDELS GOODS-S Deep Region 



ID 


RA & DEC (J2000) 


-^125 




■^125 — ^^160 


[3.6]'^ 


[4.5]'' 


Other ID'^ 


AUTO.048 


3:32:49.936 


-27:48:18.101 


25.65±0.07 


1.1±0.2 


-0.2±0.1 


C 


C 


ISO.085; 012-2499448181 


AUTO.lOO 


3:32:41.417 


-27:44:37.831 


25.98±0.14 


1.8±0.6 


0.0±0.2 


B 


B 


IS0.164 


AUT0.212 


3:32:20.965 


-27:51:37.073 


26.39±0.15 


l.OitO.3 


0.1±0.2 


> 26.8 


> 26.2 


012-2209751370 


AUT0.293 


3:32:20.981 


-27:48:53.467 


26.88±0.17 


> 1.1 


-0.4±0.3 


26.16±0.13 


26.12±0.22 


ISO.071; 012-2209848535 


AUTO_204 


3:32:18.185 


-27:52:45.566 


27.38±0.19 


0.9±0.5 


0.2±0.3 


> 26.7 


> 26.2 


012-2181852456 


AUTO.094 


3:32:40.675 


-27:45:11.624 


27.39±0.27 


0.9±0.6 


0.0±0.4 


C 


C 




AUTO.035 


3:32:34.998 


-27:49:21.623 


27.55±0.21 


> 1.2 


0.3±0.3 


26.29±0.22 


25.67±0.25 


ISO.063; 012-2350049216 



^' For the objects that arc not detected in Vios, their Vios — limits are caculated using the 2-(t flux upper limits in Yios as measured within the 

MAG_AUTQ apertures. 

^' The IRAC magnitudes arc MAG_APER within r — lV5 aperture, and limits arc based on 2-(7 flux upper limits within the same aperture, both of which 
arc with an aperture correction to a total flux. "B" or "C" means that the object is blended with or contaminated by foreground neighbor(s) in IRAC. 

For the common candidates in the AUTO and the ISO samples (see appendix), their ZD's in the other set are given. In addition, if an object is also 
in the 012 sample, its ID in 012 is given as well. 
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Fig. 6. — Similar to Figure 1, but for selection of Yios-dropouts using MG_ISO magnitudes. 



APPENDIX 

Due to the nature of the dropout selection technique, it is not unusual that samples constructed by independent 
groups can differ significantly on an object-by-object basis (see §3). To further demonstrate this point, we constructed 
a separate Kios-dropout sample to mimic an independent study. 

We used exactly the same data as in §2, and therefore any difference in this sample should be attributed to the 
selection technique. We followed the same procedures as in §3, but used MAG_ISO instead when calculating colors. The 
motivation of adopting MAG_ISO in this exercise is that isophotal apertures are commonly used in measuring Lyman- 
break galaxy colors (e.g., Steidel et al. 2003) because they often lead to higher S/N, even if they do not necessarily 
include all of the light from each galaxy. Hereafter we refer to this sample as the "ISO sample" to separate from the 
one described in the main text, which we refer to as the "AUTO sample". The initial candidates in the ISO sample 
were also visually examined by the same group of six inspectors. The visual inspection of this sample was done about 
one month later than that for the AUTO sample, and therefore even when we were examining the same object that 
is also in the AUTO sample we did not keep the memory of the result from the last time, and hence we mimicked an 
"independent" study as much as we could. 

The final ISO sample consists of eleven candidates. Fig. 6 shows their selection in the yio5>/i25-ffi60 color space, 
while Fig. 7 shows the redshift selection function. Images of these candidates are displayed in Fig. 8. Table 2 lists 
their photometric information. 

For the sake of completeness, we also compared the ISO sample to the 012's sample. Based on the photometry 
reported in 012, nine of their candidates are expected to meet our color criteria (see also §3.3.3). However, only three 
of them are included in our ISO sample (see Table 2). Among the five that are not in our sample, 012's CANDY- 
2432246169 has S/N < 5 in J125 and Hiqq and is not included in our MAG.ISD catalog. CANDY-2272447364 has 
S/N > 2 in VeoG and does not satisfy our optical non-detection criterion. CANDY-2277945141 and 2320345371 have 
Y105 — J125 — 0.74 and 0.57 mag based on our photometry and hence do not meet out color criteria. The other two, 
012's CANDY-2209751370 and 2181852456, meet all our color criteria, but were rejected during the visual inspection 
when constructing our ISO sample. Note that these latter two objects actually are in our AUTO sample (AUT0_212 
and 204, respectively) and survived the visual inspection for our AUTO sample (see below). 

A more constructive comparision is that between this ISO sample and the AUTO sample, which shows that they 
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Fig. 7. — Similar to Figure 2, but for the ISO Yios-dropout sample. The labeled magnitude ranges are still based on MAG_AUTO in J125, as 
the simulation only has control over the input total magnitudes, which are best represented by MAGJVUTQ. To be consistent, the correction 
of V^f f (m) for the ISO sample is done at the corresponding MAG_AUTO ranges. 
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Fig. 8. — Similar to Fig. 2, but for the ISO sample. The ones that are in common with the AUTO sample are not displayed here. 



have four objects in common (see Table 2). There are seven ISO candidates not in the AUTO sample. Five of them 
were rejected in the AUTO selection because they have S/N > 2 within their MAG_AUTO apertures in at least one veto 
image (IS0_011, 017, 157, 158, 160). One was rejected because its MAG_AUTD Y — J color limit does not satisfy 
F - J > 0.8 mag (ISD_008). The other one (IS0_082) satisfies ah the quantitative criteria in the AUTO selection, 
however it was rejected in the visual inspection step during the AUTO selection run. 

On the other hand, four AUTO candidates are not in the ISO sample. Two of them were rejected because they have 
MAG_AUTO Y - J color or limit below the F - J > 0.8 mag threshold (AUT0_094 and 368). The other two satisfy ah 
the quantitative criteria in the ISO selection, however they were rejected in the visual inspection step during the ISO 
selection run (AUT0_204 and 212). 

This internal comparison of our two samples thus further demonstrates the points addressed earlier: (1) adoption of 
different photometry result in samples that can be significantly different on an object-by-object basis, and (2) visual 
inspections at different time, especially the inspections of the veto images, can also result in differences in this sense 
because such inspections are not guaranteed to be fully repeatable for a large sample at the S/N< 2 level. 

Nevertheless, the main statistical trends revealed by the ISO sample are very similar to those inferred from the 
AUTO sample. Both samples include the same two bright Fio5-dropouts (AUT0_048/IS0_085 and AUTD_100/IS0_184), 
and the same two IRAC-detected ones (AUT0_293/IS0_071 and AUTD_035/ISD_063). For completeness. Fig. 9 shows 
the constraint from the ISO sample on the LF, where one can see that the bright-end excess is still present. Similarly, 
the constraints on the MF and the stellar mass density derived from the ISO sample agree with those based the AUTO 
sample to within < 10%, and the small differences are mainly caused by the slightly different Ug// corrections over the 
range of interest. From this exercise, we believe that the main results presents in this work are robust. 
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Fig. 9. — Same as Figure 5, but for the ISO sample. Note that the J125 magnitudes shown here are the MAGJVUTO magnitudes, because 
these are taken as the total magnitudes. 



TABLE 2 

MAG.ISD-SELECTED Y105-DROPODTS IN CANDELS GOODS-S Deep Region 



ID 


RA & DEC {J2000) 




I105 - 


J125 - Hiao'' 


[3.6]d 


[4.5]d 


Other ID<= 


ISO.085 


3:32:49.936 -27:48:18.100 


25.65±0.07 


1.2±0.1 


-O.litO.l 


C 


C 


AUTO.048; 012-2499448181 


IS0.164 


3:32:41.417 -27:44:37.831 


25.98±0.14 


0.9±0.1 


O.OitO.l 


c 


c 


AUTCIOO 


IS0.157 


3:32:42.882 -27:45:04.268 


26.50±0.13 


0.9±0.2 


0.1±0.2 


B 


B 




ISO.071 


3:32:20.981 -27:48:53.468 


26.88±0.17 


1.7±0.8 


0.1±0.2 


26.16±0.13 


26.12±0.22 


AUT0.293; 012-2209848535 


ISO.082 


3:32:14.133 -27:48:28.911 


27.16±0.19 


1.1±0.5 


0.0±0.2 


> 27.1 


> 26.5 




ISO.Oll 


3:32:14.469 -27:51:48.542 


27.24±0.30 


1.2±0.6 


0.0±0.3 


B 


B 




IS0.158 


3:32:47.953 -27:44:50.436 


27.28±0.23 


0.8±0.3 


0.2±0.2 


> 26.7 


> 25.9 




ISO.063 


3:32:34.999 -27:49:21.622 


27.55±0.21 


> 1.5 


0.1±0.2 


26.29±0.22 


25.67±0.25 


AUTO_035; 012-2350049216 


ISO.017 


3:32:18.091 -27:51:18.492 


27.57±0.18 


1.3±0.7 


-0.4±0.3 


B 


B 




ISO.160 


3:32:46.111 -27:44:47.997 


27.88±0.50 


0.9±0.4 


0.1±0.3 


> 26.7 


> 26.1 




ISO.008 


3:32:16.915 -27:52:01.878 


28.00±0.27 


> 0.9 


-0.3±0.4 


> 26.7 


> 26.2 





^' The quoted total magnitudes in J125 are their MAG_AUTO values. 

^' The Yio5 — J125 colors are based on their MAG.ISO magnitudes. For the objeets that are not detected in Vios, their Yios — J125 limits are eaculated using 
the 2-{T flux upper limits in Yios as measured within the MAG.ISO apertures. 

The J125 — ^iQQ colors arc based on their MAG_ISO magnitudes. 
^' The IRAC magnitudes are MAG_APER within r — l''5 aperture, and limits are based on 2-(T flux upper limits within the same aperture, both of which are 
with an aperture correction to a total flux. "B" or "C" means that the object is blended with or contaminated by foreground neighbor{s) in IRAC. 

For the common candidates in the AUTO and the ISO samples, their ID's in the other set are given. In addition, if an object is also in the 012 sample, 
its ID in 012 is given as well. 



